Orbital ordering in e g orbital systems: Ground states and thermodynamics of the 120° model 
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Orbital degrees of freedom shape many of the properties of a wide class of Mott insulating, transition metal 
oxides with partially filled 3<i-shells. Here we study orbital ordering transitions in systems where a single 
electron occupies the e 9 orbital doublet and the spatially highly anisotropic orbital interactions can be captured 
by an orbital-only model, often called the 120° model. Our analysis of both the classical and quantum limits of 
this model in an extended parameter space shows that the 120° model is in close proximity to several T = 
phase transitions and various competing ordered phases. We characterize the orbital order of these nearby phases 
and their associated thermal phase transitions by extensive numerical simulations and perturbative arguments. 



PACS numbers: 71.20.Be, 05.70.Fh, 75.25.Dk 

Mott insulating transition metal oxides with partially filled 
3e?-shells - such as the manganites - exhibit rich phase dia- 
grams with many competing orders, indicating a non-trivial 
interplay of spin, charge, and orbital degrees of freedom JT]. 
A prominent example of a material exhibiting orbital order is 
the extensively studied LaMnC>3 10. The crystal field in this 
perovskite material splits the five rf-orbitals into three l 2 g or- 
bitals occupied by three electrons, and an e g doublet sharing a 
single electron. This partially filled e g doublet then gives rise 
to an additional orbital degree of freedom indicating which of 
the two orbitals is occupied. The exchanges between these 
orbital degrees of freedom - arising from Jahn-Teller distor- 
tions or Kugel-Khomskii type superexchange - are oftentimes 
described by orbital-only models which neglect the spin de- 
grees of freedom. The latter is justified in situations where the 
energy scales of spin and orbital interactions are well sepa- 
rated, i.e. orbital interactions correspond to temperature scales 
where the spins are still largely disordered or in a situation 
where the spins are effectively frozen out (e.g. by a magnetic 
field). Expressing the e g orbital degree of freedom by a two- 
component pseudospin T = (T z ,T X ), where T z — ±1 corre- 
spond to occupation of the |3z 2 — r 2 ) and \x 2 — y 2 ) orbitals, 
the highly anisotropic interactions between them are captured 
by the so-called 120° model (4) on a cubic lattice 

#120 = - J2i n =x,y \ \j^ T l T Ul + : '>-l, / ,' / ;'. - (1) 

±V3J mix (77 ?7 +7 + Tf Tf +T )] - J Z T Z T Z +Z , 

where the ± sign for the 'mixing' term enters for coupling 
along the x and y directions, respectively Q. If the or- 
bital exchange is primarily mediated through Jahn-Teller dis- 
tortions, this model is commonly considered in its classi- 
cal limit, where the pseudospins T are 0(2) spins. If, on 
the other hand, the orbital exchange arises primarily from a 
Kugel-Khomskii type superexchange 0, this model should 
be considered in its quantum limit. In the latter case, the 
pseudospins T are identified with SU(2) spins, i.e. their 
components become Pauli matrices T x,z — \o x,z ■ The 
above 120° model has typically been studied at equal cou- 



pling J x = J z = </ m ix, for which it exhibits an enhanced 
rotational symmetry where the symmetry of the cubic lattice 
under permutations of the x, y, and z axes is reflected in a 
three-fold symmetry in the (T Z ,T X ) plane. This becomes 
apparent when rewriting {[} as H 120 = -J ~Y,i,~ t = x ,y,z( T >- ' 
e 7 )(ri+ 7 • e 7 ) , where the e 7 are unit vectors in the x, y, z- 
directions and the Tj are defined as three-component vectors 
n = {[T( + V3Zf]/2, [Tf - VZT?]/2, Tf). 

While the presence of this enhanced rotational symme- 
try for equal coupling has greatly benefitted the understand- 
ing of the classical model and has led to a rigorous descrip- 
tion of its highly degenerate ground-state manifold |6|, it has 
remained elusive to identify the ground states of the quan- 
tum model solely based on symmetry arguments. In this 
manuscript, we will take a broader perspective and study the 
above 120° model away from this symmetric point and ex- 
plore ground states and thermodynamic properties in an ex- 
tended two-dimensional parameter space (J x / J z , J m \^/ J z ), 
which in experiments should be accessible by changing pres- 
sure or adding a small electric field. Our approach reveals 
that the original 120° model is in close proximity to sev- 
eral T = phase transitions and various competing ordered 
phases. Combining extensive numerical simulations with ana- 
lytical arguments we describe the orbital order in these phases 
for the classical and quantum limits of this extended 120° 
model, as well as thermal phase transitions associated with 
these phases and T — phase transitions between them. 

The classical model- For the classical 120° model with 
rotational symmetry, e.g. J x = J z = J m i x , it has long 
been appreciated that this model exhibits an infinite, but sub- 
extensive ground-state degeneracy [6], which is split at low 
temperatures by an order-by-disorder mechanism stabilizing 
six ordered states |6]|7]- Before turning to the question of 
how these characteristic features change when exploring the 
model in the extended parameter space, we will briefly re- 
count their origin in the symmetric model. To this end, we 
label a general state in the orbital subspace by an angle \0) = 
cos(0/2) \3z 2 - r 2 )+sm(f?/2) \x 2 - y 2 ), which for the clas- 
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sical model simply describes the orientation of the 0(2) pseu- 
dospin vector. To identify the degenerate manifold of ground 
states, we first observe that any polarized state with all pseu- 
dospins being aligned along some angle 9* is a ground state of 
Hamiltonian ([TJ. Starting from any such state, further ground 
states can be found [6] [7] by reflecting all orbitals in the xy 
plane about a line at 0°, the xz plane about 120°, or the yz 
plane about 240°. Remarkably, we find that this ground-state 
degeneracy remains (partially) unscathed for an extended pa- 
rameter regime when moving away from the symmetric model 
on a line described by < J m ix/Jz < 1 and J x — J z . What 
distinguishes states along this line in parameter space, how- 
ever, is their instability to thermal fluctuations and the entropic 
selection of low-temperature states. To discuss this order- 
by-disorder mechanism we calculate the free energy of the 
low-temperature states by considering a spin-wave approxi- 
mation of ([TJ and expanding to second order in small fluc- 
tuations 89i = 6i — 9* about an orbitally ordered state with 
9i = 9* at each site. The resulting free energy F(9*) is plot- 
ted in Fig. [T] as a function of J m i X /J z . For J mix < 0.8 J z , 
the ground-state manifold is lifted and four low-temperature 
states are entropically favored with their free energy being 
minimized at angles 9* = 0°, 90°, 180°, 270°. These four 
states correspond to orbitally ordered states in orbital configu- 



rations given by |0°) = 3z 



2 ), |180° 
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y 2 )) /y/2 for the 90° and 270° states, 



respectively. For J m j x > 0.8 J z the 90°-minima in the free- 
energy curves bifurcate and in total form six minima, all of 
which become exactly equal only for the symmetric model 
Jmix = Jz, where the minima are located precisely at angles 
of 9* = 0°, 60°, 120°, . . . , 300°. If we further enlarge J mix 
beyond J z , this order-by-disorder phenomenon disappears 
and we instead find that the ground and low-temperature states 
are energetically selected, with the xz or yz planes ordering 
in alternating orientations of 9\ and 9\ + 180° where 9\ con- 
tinuously changes from 9\ w 30° to 9\ 45° with increasing 
J m ; x . This transition between entropic and energetic selec- 
tion occurs exactly at the symmetric point J x = J z = J m ; x 
of the 120° model and is accompanied by a first-order phase 
transition at zero temperature, which becomes apparent in a 
level-crossing of ground-state energies shown in Fig. [3] 
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FIG. 1 : Entropic selection of low-temperature states in the classical 
model: The free energy F{6*) obtained from a spin- wave analysis 
of {TJ as a function of J m i x / J z for fixed J x = J z . 
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FIG. 2: Orbital ordering transition in the classical model: Specific 
heat C V (T) (upper panel) and order parameter M(T) (lower panel) 
versus temperature for various system sizes L. 



We now turn to a discussion of the thermodynamic prop- 
erties of the classical 120° model, in particular the thermal 
ordering transition into the low-temperature states described 
above. To investigate the latter we have run extensive Monte 
Carlo simulations of model ([TJ, going well beyond previous 
numerics for the (diluted) symmetric model [8]. Concentrat- 
ing on a family of models, where we vary J m i x /J z but keep 
Jx = Jz fixed, we find a line of continuous thermal phase 
transitions. Fig. [2] shows the specific heat C V (T) at the tran- 
sition diverging with linear system size L for two members in 
this family, the symmetric model with J m j X = J z and a 'trun- 
cated' model where we drop the mixing terms in Hamiltonian 
([TJ, i.e. J m j x = 0. For all models in this family we can cap- 
ture the transition to the low-temperature ordered states by a 
single order parameter M (independent of J m ; x < J z ). Since 
in our numerical simulations we do not know a priori which 
one of the three possible ordering planes the system sponta- 
neously selects at the ordering transition, we define the order 
parameter as the maximum of the xy, xz, and yz plane mag- 
netizations M = max(M xy , M xz , M yz ), where the magneti- 
zation in the xy plane is given by M xy = J2 Z J2 X , V T x, v 
and M xz , M yz are obtained by cyclic permutations of the in- 
dices in the sums. As expected this order parameter quickly 
grows at the transition temperature T c (see the lower panel of 
Fig.|2j. Despite the relatively large system sizes studied here, 
finite-size effects still render the identification of the univer- 
sality class of these transitions somewhat ambiguous, remi- 
niscent of studies of similar models in two spatial dimensions 
||9l . Tracking the ordering temperature T c with the strength of 
the mixing term J m i X , as shown in Fig. [3] we find a signifi- 
cant suppression for the symmetric 120° model, for which the 
transition occurs around T c /J z — 0.677 ± 0.003. 

The quantum model- As an inroad into exploring ground 
states and thermodynamics of the quantum 120° model in an 
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FIG. 3: Variation of the transition temperature (upper panel) and 
ground-state energy (lower panel) in the classical model for varying 
Jmix- The 120° model corresponds to J m i x = J z . 



extended (J x /J z , Jmix/J z ) parameter space around the sym- 
metric point, we first concentrate on a family of models for 
which the mixing term J mix vanishes but for which we can 
still vary J x /J z . This line in parameter space stands out 
as it allows for a thorough analysis using quantum Monte 
Carlo (QMC) simulations, while all other regions in param- 
eter space of non-zero J m i x are plagued by the so-called sign 
problem. Along this line, we first discuss the 'truncated' 
model at equal coupling J. x = J z , for which we have run 
extensive QMC simulations using an extension of the ALPS 
looper code iflOl fTTI . Our numerical findings, summarized 
in Fig. [4] show that this model undergoes a continuous ther- 
mal phase transition around T c /J z = 0.41 ±0.01 into an 
orbitally ordered state at low temperatures. In this ordered 
state all orbitals are found to spontaneously orient in either the 
|3z 2 — r 2 ) or \x 2 — y 2 ) orbital configurations, correspond- 
ing to pseudospins pointing in the ±T 2 directions as indi- 
cated in Fig. |4j5). This ordered orbital state, which we call the 
'T z polarized' state, precisely corresponds to the 0° and 180° 
states found as low-temperature states in the classical trun- 
cated model, indicating that thermal fluctuations and quantum 
effects favor the same states. 

As we vary J x /J z away from equal coupling we find that 
the truncated model above exhibits another peculiarity: it 
sits right at a first-order transition between different quan- 
tum ground states ifT^l . For J x < J z this is the same 
±T Z polarized state found for equal coupling (at finite tem- 
perature), while for J x > J z the pseudospins in any given 
xy-plane point along the ±T X directions corresponding to 



(|3z 2 -r 2 ) ± 



V 



! ))/-\/2 orbital configurations (illus- 



trated in the insets of Figs. [5] and [6), but pseudospins in dif- 
ferent xy-planes do not have to be aligned. This first-order 
transition is apparent in the level crossing of energies shown 
in the lower panel of Fig. [5]calculated from both QMC simu- 
lations (at temperatures well below the thermal transition) and 
2nd order T = perturbation expansions around the limits of 
J x = and J x — > oo. The rather good agreement of the 
two approaches indicates that quantum effects only modestly 
change the ground states with varying J x /J z . 



Both orbitally ordered phases discussed above exhibit 
gapped elementary excitations corresponding to a single pseu- 
dospin flip, e.g. an 'orbital flip' |3z 2 — r 2 ) <H- \x 2 — y 2 ) in 
the T z polarized state. We can directly estimate the excita- 
tion gap A(T) of such an 'orbiton' excitation in our QMC 



>P simulations (TT] El as A(T) = 2nT 
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where S(uj) is a Fourier transform on an imaginary time cor- 
relation function ^ r C(r, r) of the pseudospins. Our results 
are given in the top panel of Fig. [6] Again we find good quan- 
titative agreement with 3rd order perturbative results (dashed 
lines) calculated around the limits of J x = 0,oo. For both 
orbitally ordered phases, the orbiton gap is suppressed as one 
approaches the truncated model at equal coupling, where for 
both phases we measure a gap of A( J x — J z ) = (0.34 ± 
0.04) J z . As a consequence, the temperature of the thermal 
phase transition is also significantly suppressed in the vicinity 
of the truncated model (see Fig. [5]). 

Having established the nature of the ordered phases for van- 
ishing mixing term, we are now in a position to return to a 
discussion of the 120° model in the full ( J x j J z , J m i X / J z ) pa- 
rameter space. While a small mixing term J m j X =^ will 
not affect these gapped phases, a sufficiently large mixing 
term can close the orbiton gap. For large J m ; x , a mean-field 
approximation indicates an 'orthogonal' ordered orbital state 
where the pseudospins are aligned within each xz (or yz) 
plane, but aligned approximately perpendicular to each other 
between planes. Assuming that the instability of the polar- 
ized phases indeed arises primarily from orbiton condensa- 
tion (and is not preempted by some other transition), we can 
map out a phase boundary in parameter space by calculating a 
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FIG. 4: Orbital ordering transition in the quantum model: a) specific 
heat C V (T), b) magnetization M(T) = j £\ T*\, and c) Binder cu- 
mulant B [M(T)\ = 1 - <A/(T) 4 ) /3(M(T) 2 } 2 . The latter exhibits 
a crossing point for data of different system sizes L, strongly indica- 
tive of a continuous transition [ 13] • 
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FIG. 5: Variation of the transition temperature (upper panel) and 
ground-state energy (lower panel) for the quantum model. 



T = perturbation expansion of the orbiton gap in J m i x /J z . 
Our results are shown in the lower panel of Fig. [6| where the 
lines indicate the closing of the orbiton gap when consider- 
ing 2nd order corrections to the gap values for a given ratio 
J x /J z calculated either from QMC simulations (symbols) or 
perturbation theory (dashed lines). Interestingly, these phase 
boundaries intersect the J x = J z axis at values of J m i X /J z 
just above one, the location of the symmetric 120° model. 
Given that the perturbative results overestimate the critical 
value of J m i X /J z , one might be tempted to conclude that 
the symmetric model is right at a multicritical point between 
the three phases in this extended phase diagram. In particu- 
lar, this would indicate that the quantum ground states found 
for the truncated model at J m [ x /J z = adiabatically con- 
nect to the ground states at the symmetric point. However, at 
the symmetric point the rotational symmetry requires that if 
the |3z 2 — r 2 ) (|x 2 — y 2 )) orbital states (i.e. the ±T Z polar- 
ized states) remain ground states, then also their symmetry re- 
lated \3x 2 — r 2 ) (\y 2 — z 2 )) and |3y 2 — r 2 ) (|z 2 — x 2 )) or- 
bital states must become ground states. This leaves us with 
two possible scenarios to connect the T x polarized states of 
the truncated model to the symmetric point: i) Quantum ef- 
fects for non-zero J m j x have the same effect as thermal fluc- 
tuations in the classical model and these states adiabatically 
turn into a combination of the symmetry required states above 
(with two more states coming down from higher energies). 
Some support for this scenario comes from considering a 1/S 
expansion for varying < Jmix/Jz < 1 (generalizing previ- 
ous calculations [ 14, 15 1), which in linear order gives a zero- 
point energy that exactly mimics the behavior of the free en- 
ergy obtained for the classical model in Fig.[T] i.e. two minima 
at 90° and 270° found in the vicinity of the truncated model 
bifurcate with increasing J m i x /J z . ii) The T x polarized state 
remains unchanged and brings in four more symmetry related 
states. Scenario i) thus selects exactly the same six order- 
ing states at angles 0° , 60° , ... as in the classical symmetric 
model, while for scenario ii) we get twelve ordering states at 
angles 0°, 30°, 60° ...,330°. 

We acknowledge discussions with L. Balents, A. Laeuchli, 
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FIG. 6: Phase diagram of the quantum 120° model: The upper panel 
shows the orbiton gap of the orbital polarized states as a function 
of Jx/Jz obtained from QMC simulations at T = 0.15 (symbols) 
and 3rd order perturbation theory (dashed lines). The lower panel 
shows the critical value of J m j x at which the orbiton gap closes. The 
location of the symmetric 120° model is indicated by the arrows. 
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